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Human papillomavirus type 16 (HPV16) and other oncogenic viruses have been reported to 
deregulate immunity by suppressing the function of the double-stranded DNA innate sensor 
TLR9. However, the mechanisms leading to these events remain to be elucidated. We show 
that infection of human epithelial cells with HPV1 6 promotes the formation of an inhibi- 
tory transcriptional complex containing NF-KBp50-pG5 and ERa induced by the E7 onco- 
protein. The E7-mediated transcriptional complex also recruited the histone demethylase 
JARID1B and histone deacetylase HDAC1. The entire complex bound to a specific region on 
the TLR9 promoter, which resulted in decreased methylation and acetylation of histones 
upstream of the TLR9 transcriptional start site. The involvement of NF-kB and ERa in the 
TLR9 down-regulation by HPV1G E7 was fully confirmed in cervical tissues from human 
patients. Importantly, we present evidence that the HPV1 6-induced TLR9 down-regulation 
affects the interferon response which negatively regulates viral infection. Our studies 
highlight a novel HPV1 6-mediated mechanism that combines epigenetic and transcriptional 
events to suppress a key innate immune sensor. 



Approximately 26% of all human cancers world- 
wide are associated with infectious agents, among 
which 80% are viruses (Bouvard et al., 2009). 
Hepatitis C virus (HCV), hepatitis B virus (HBV), 
Human T-lymphotropic virus type I (HTLV-1), 
Epstein-Barr virus (EBV), Kaposi sarcoma- 
associated virus (KSHV), and the mucosal high- 
risk (HR) human papillomavirus (HPV) types 
have been clearly implicated in different types 
of human cancers. In particular, HPV is the eti- 
ological factor of cervical cancer and is respon- 
sible for ^20% of all human cancers linked to 
infection. More than 500,000 new cervical can- 
cer cases and 275,000 deaths are reported each 
year worldwide (Tay, 2012). In addition, HR 
HPV types are responsible for a proportion of 



oropharyngeal cancers (Ryerson et al., 2008; 
Marur et al., 2010) that appear to have steadily in- 
creased in the last two decades in the USA and 
Europe (Ryerson et al, 2008; Marur et al, 2010). 
Similarly to several oncogenic viruses, HR HPV 
types are able to alter the immune surveillance 
and cellular homeostasis, mainly deregulating the 
cellular gene expression and promoting epigen- 
etic changes (Kim et al., 2003; Zhang et al., 2004; 
Rincon-Orozco et al., 2009; Moore and Chang, 
2010; Nelson, 2011). The initial outcome is the 
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increase of cellular proliferation that favors viral replication 
and persistence. This event is mainly mediated by the HR 
HPV oncoproteins E6 and E7, which have the properties 
to target many cellular proteins. In particular, E6 and E7 
have the properties to bind and degrade the tumor suppressor 
gene products p53 and pRb, respectively. As a consequence, 
the transcription of several cell cycle check point tumor sup- 
pressors, such as p21, FAS, DEK, and B-MyB, is prevented 
(Tommasino et al., 2003; Malanchi et al., 2004; Rampias et al., 
2009). In addition, the HR HPV types are able to deregulate 
several immune-related pathways to guarantee the persistence 
of the infection. 

Innate immunity is the first line of host defense against in- 
fections. The cellular innate immune response is mediated by 
both hematopoietic and nonhematopoietic cells that express 
pathogen recognition receptors (PRRs). Viral nucleic acids 
are sensed by different PRR families, which include: cyto- 
solic retinoic acid-inducible gene I (RIG-I)— hke helicases for 
viral RNA; cytosolic DNA sensors such as DNA-dependent 
activation of IFN (DAI), Absent in Melanoma 2 (AIM2), and 
others; and several members of the TLR family (Coban et al., 
2005; Adams et al, 2008; Fernandes-Alnemn et al, 2009; 
Gondois-Rey et al., 2009; Rathmam et al, 2010). TLR9 was 
one of the first innate immune receptors to be identified 
(Hemmi et al., 2000). Historically, in humans TLR9 was re- 
ported to be expressed in plasmacytoid DCs (pDCs) and B cells. 
However, recent studies have shown that it is also present in 
nonimmune cells such as endothelial and epithelial cells (Lebre 
et al., 2007; Pivarcsi et al, 2007; Monzane et al., 2012). Upon 
recognition of unmethylated double-stranded DNA CpG 
motifs present in the genome of viruses such as EBV, HSV, 
and HPV (Lund et al, 2003; Hasan et al, 2007a; Fiola et al, 
2010; Zauner et al., 2010), TLR9 initiates a signaling cascade 
that leads to the production of type I IFN and proinflamma- 
tory cytokines (Sepulveda et al., 2009; Sasai et al., 2010; Avalos 
and Ploegh, 2011; Ewald et al., 2011), the release of which 
activates host immune defenses against infection. 

Despite the efficiency of innate immune response, HR 
HPV types have developed strategies to persist in the host. 
Indeed, HR HPV16 and 18 are able to inhibit the transcrip- 
tion of proinflammatory chemokines and cytokines, such as 
CCL5 and IL-1(3 (Karim et al, 2011). In addition, the IFN 
inducible antiviral genes IFIT1 and MX1, proapoptotic genes 
(TRAIL and XAF1), and PRRs (TLR3, RIG-I, and MDA5) 
are also inhibited by HPV16, 18, and 31 (Reiser et al, 2011). 
We and others have observed that HPV16, EBV, and HBV 
impair the expression and function of the innate immune re- 
ceptor TLR9 (Hasan et al, 2007a; Fathallah et al, 2010; Hirsch 
et al., 2010;Vincent et al., 2011). TLR9 expression is severely 
suppressed in the cervical epithelium of women with HPV16- 
positive cervical lesions compared with that of women that 
are healthy or with regressing infection (Hasan et al., 2007a; 
Daud et al., 2011). The viral oncoproteins E6 and E7 have 
been linked to the HPV16-mediated TLR9 transcriptional 
down-regulation (Hasan et al., 2007a). However, the mecha- 
nisms involved in this event remain to be elucidated. 



Here, we have elucidated the mechanism of HR HPV16 
E7 mediated down-regulation of TLR9. This viral oncoprotein 
induces the formation of transcriptional inhibitory complexes, 
including NF-KBp50— p65, ERa, and chromatin modifying 
enzymes, at the TLR9 promoter and induce epigenetic changes. 
Thus, HPV1 6 employs a unique mechanism to turn downTLR9 
transcription, expression, and function, which is an essential 
event required for oncoviral mediated carcinogenesis. 

RESULTS 

Infection of human epithelial cells with HPV16 
down-regulates TLR9 expression and function 
in an E6/E7-dependent manner 

To understand how TLR9 expression is regulated by HPV16 
we generated quasi-virions (16QsV) that closely resemble the 
natural virus (Flores et al., 1999; Pyeon et al., 2005). Infection 
of C33A cells (an HPV-negative cervical epithelial cell line) 
with 16QsV for 8 h led to reduced TLR9 mRNA levels with 
a further decrease observed at 24 h (Fig. 1 A). The effect seen 
at 24 h after infection was dependent on viral genome repli- 
cation, as UV-treated 16QsV did not suppress TLR9 levels 
(Fig. 1 A, right). In contrast, CpG 2006 and HSV-2, which 
strongly activate the TLR9— NF-kB axis, transiently suppressed 
TLR9, with restoration of mRNA and protein levels 24 h 
after stimulation (Fig. 1 A). Addition of TNF, pseudo-virions 
(late proteins LI and L2 encapsidated GFP expression plas- 
mid, abbreviated as PV), and GpC controls did not down- 
regulate TLR9 levels (Fig. 1 A). Thus, TLR9 expression was 
severely suppressed by infection with 16QsV. 16QsV infec- 
tion and replication was controlled by qPCR using specific 
primers for the early genes El and E7 (Fig. 1 B, left and 
middle). The number of added 16QsV viral genome equiva- 
lents inversely correlated with TLR9 mRNA levels (Fig. 1 B, 
right). We previously reported that TLR9 down-regulation 
was associated with E6 and E7 expression in primary kerati- 
nocytes as well as in cervical cancer— derived cell lines (Hasan 
et al., 2007a). Accordingly, C33A cells infected with 16QsV for 
24 h in the presence of a siRNA against HPV16E6E7 restored 
TLR9 mRNA levels and promoter activity (Fig. 1 , C and D) . 
We next investigated the biological consequence of TLR9 
suppression by HPV16. Human keratinocytes are strong pro- 
ducers of proinflammatory cytokines such as IL-6, IL-8, and 
MIP3a (Debenedictis et al., 2001; Hudak et al., 2002; Ito et al, 
2003; Metz et al., 2008; Bangert et al, 2011; Kaplan et al., 
2012). We have previously shown that HPV16 E6E7 prevented 
secretion of MIP3a and IL-8 when cells were stimulated as a 
result of the loss of TLR 9 expression (Hasan et al., 2007a). 
We next evaluated whether infection with 16QsV also sup- 
pressed TLR9 functional signaling. C33A cells were infected 
with 16QsV for 36 h, washed, and stimulated with a CpG 
motif from HPV 16 genome (Hasan et al, 2007a) or CpG 2006. 
We observed that 16QsV infection hindered TLR9 func- 
tion, as CpG 2006 and CpG motifs from HPV16 did not lead 
to the secretion of IL-8, IL-6, and MIP3a (unpublished data). 
TLR9 is also a strong inducer of type I IFN, the release of 
which activates host immune defenses against viral spread 
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Figure 1 . TLR9 expression is suppressed 24 h after infection by native 1 6QsV. (A, left) C33A cells were not treated (NT) or treated for 8 or 24 h 
with TNF, PV, 16QsV (at 10 7 viral concentrations genome equivalents), 16UV (rendered replication incompetent using UV), CpG 2006, GpC 2006, or in- 
fected with HSV-2. TLR9 mRNA levels were determined by qPCR. Shown are data from four independent experiments performed in triplicate. Error bars 
indicate SEM. (A, right) TLR9 protein was examined by immunoblotting in C33A cells. Cells were harvested after 24 h treatment with PV, TNF, 16UV, 16QsV 
CpG 2006, and GpC control. (B) C33A cells were treated with increasing viral concentrations genome equivalents (v.g.e.; as measured by qPCR on the viral 
DNA of infected cells) of 16QsVfor 8 or 24 h). E1, E7 (left and middle), and TLR9 mRNA levels (right) were determined for their relative expression by 
qPCR. Shown are data from five independent experiments performed in duplicate with ***, P < 0.0001, based on an unpaired Student's ttest. (C) C33A 
cells were infected as indicated for 24 h. siRNA against HPV16E6E7 was transfected 24 h after stimulation and TLR9 mRNA levels were determined by 
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(Sepulveda et al., 2009; Sasai et al., 2010; Avalos and Ploegh, 
2011; Ewald et al, 2011). We tested the ability of HK trans- 
duced with HPV16E6E7 or with empty vector (PLXSN) to 
produce type I IFN upon TLR9 stimulation with CpG 2216. 
HPV16E6E7 expression severely impaired the ability of TLR9 
to produce type I IFN compared with cells transduced with 
the vector alone (Fig. 1 E, left). The same block in IFN pro- 
duction was observed in C33A cells infected with 16QsV 
before CpG stimulation (Fig. 1 E, middle) and correlated to 
the loss TLR9 mRNA levels. No effect on IFN production 
was observed when PV was used as a control Fig. 1 E, mid- 
dle). Addition of CpG 2216 24 h before 16QsV infection 
permitted type I IFN activation of the ISRE minimal promoter 
that was abrogated in the presence of an antibody against the 
type I IFN receptor (anti-IFNR; Fig. 1 E, right). Most im- 
portantly, prestimulation of TLR9 with CpG 2216 before 
16QsV infection significantly decreased the expression of the 
viral early genes which was blocked in the presence of anti- 
IFNR (Fig. 1 F). To demonstrate that this event was linked 
to TLR9 down-regulation and not to the alteration of IFN 
signaling (Ronco et al., 1998), we tested whether 16QsV blocks 
the type I IFN production signaling pathway of RIG-I. Indeed, 
ectopic levels of RIG-I— stimulated cells infected with 16QsV 
did not affect type I IFN bioactivity (unpublished data). To gain 
more insights on the biological significance of HPV-induced 
TLR9 down-regulation, we silenced its expression in HK by 
using a short hairpin RNA (Fig. 1 G, left). Subsequently, these 
cells were infected with 16QsV and viral load was determined. 
Cells expressing TLR9 shRNA had a higher copy number of 
HPV16 genome in comparison with mock cells (Fig. 1 G, mid- 
dle). Accordingly, HPV16 viral transcription was increased in 
cells with reduced TLR9 expression (Fig. 1 G, right). Collec- 
tively, these data show that infection with 16QsV of human 
epithelial cells inhibited TLR9 expression and signaling in an 
E6- and/or E7-dependent manner and thatTLR9 plays an 
essential role in limiting HPV16 life cycle. 

HPV16 down-regulation is dependent on NF-kB signaling 

NF-kB signaling was shown to regulate TLR9 (Takeshita et al., 
2004) and we reported that deletion of putative NF-kB sites 



in the TLR9 promoter restored its transcriptional activity in 
the presence of HPV16 E6 and E7 (Hasan et al., 2007a). We 
next determined whether TLR9 down-regulation induced by 
16QsV is mediated by the NF-kB pathway. C33A were tran- 
siently transfected with the TLR9 promoter/luciferase reporter 
gene and treated with siRNA for IKKa or IKK(3 (Fig. 2 A, 
right), two cytoplasmic kinases that promote the nuclear trans- 
location of active NF-kB transcriptional factor (Hacker and 
Karin, 2006). Cells were then exposed for 24 h to 16QsV or 
TNF. In the presence of IKKa or IKK(3 siRNA, TLR9 pro- 
moter activity and mRNA levels were rescued compared with 
cells treated with scramble siRNA (Fig. 2 A). Interestingly, 
TNF, a known activator of the NF-kB pathway, was unable 
to inhibit TLR9 transcription (Fig. 2 A). Ectopic expression 
of a dominant-negative MyD88 mutant did not restore TLR9 
transcription or protein levels in cells infected with native 
16QsV, indicating that a MyD88-NF-KB pathway was not 
involved in this phenomenon (unpublished data) . In contrast, 
the suppression of TLR9 expression by UV-treated 16QsV 
or HSV2, which both contain CpG elements (Hasan et al., 
2007a), was prevented in the presence of a dominant-negative 
MyD88 mutant. 

A 1-h treatment with a chemical inhibitor of IKKB (Bay 1 1) 
also restored TLR9 mRNA and protein levels in all cervical 
cancer— derived cell hnes (Fig. 2 B).TLR9 expression upon 
Bay 11 treatment correlated with loss of NF-KBp65 nuclear 
localization (Fig. 2 C). In addition, gene silencing of IKKa, 
IKK(3, or NF-KBp65 in SiHa cells by siRNA resulted in the 
recovery of TLR9 expression, as measured by luciferase activ- 
ity or by the endogenous TLR9 levels (unpublished data). Thus, 
TLR9 transcriptional inhibition depends on the activation of 
NF-kB signaling after infection with 16QsV. 

We next characterized which HPV16 oncoprotein was 
responsible for NF-KB-dependent-TLR9 down-regulation. 
Human primary keratinocytes (HK) were transduced with 
E6 and/ or E7 and the pLXSN vector control. Immunoblot- 
ting showed that several positive regulators of the canonical 
NF-kB signaling, i.e., IKKB, p50, and p65, were activated by 
E7, and to a lesser extent by E6. Stimulation of the canoni- 
cal NF-kB pathway leads to IKK complex activation, which 



qPCR. Shown are data from four independent experiments performed in triplicate and error bars indicate SEM with ***, P < 0.0001, based on an 
unpaired Student's t test. (D) C33A cells were transfected with the TLR9 luciferase promoter and, 24 h later, treated with 16 or control PV. siRNA against 
HPV16E6E7 was transfected 24 h after infection and TLR9 luciferase levels were measured 18 h later. Shown are data from four independent experiments 
performed in triplicate, and error bars indicate SEM with ***, P < 0.0001, based on an unpaired Student's f test. (E, left) HK transduced with the pLXSN 
vector control or HPV16E7 were stimulated for 24 h with TLR9 agonist CpG 0DN 2216 (type A IFN inducer) and the IFN bioassay was performed. (E, mid- 
dle) C33A cells were unstimulated or infected with 16QsV or PV for 36 h, and cells were washed and then treated with CpG 2216 for 24 h. Supernatants 
were collected and the type I IFN bioassay was performed. (E, right) C33A cells were stimulated with GpC or CpG 2216 for 24 h and then infected with 
16QsV for 36 h. Supernatants were harvested and treated with anti-IFNR or IgG control and the type I IFN bioassay was performed as described. (F) C33A 
cells were unstimulated and infected with 16QsV for 36 h + CpG 2216 stimulation for 24 h or vice versa ± anti-IFNR or IgG control. Cells were harvested 
for RNA extraction and qRT-PCR was performed for E1 (left) or E7 (right) relative gene expression. Shown are data from four independent experiments 
performed in triplicate, and error bars indicate SEM with ***, P < 0.0001, based on an unpaired Student's r test. (G) HK cells were transfected for a total of 
48 h with shTLR9 and shControl (scramble sequence). 20 h after transfection, HK cells were infected for 24 h with 16QsV at a viral load of 10 7 v.g.e. Viral 
load was measured by qPCR using E6- and E7-specific primers. HPV16 E6 and E7 transcripts were also measured by qPCR and TLR9 transcripts were 
assessed by qPCR after 24 h transfection with shRNA. Data are representative of six independent experiments performed in triplicate. Shown are the mean ± 
SEM with ***, P < 0.0001, based on an unpaired Student's f test. 
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Figure 2. HPV1 6E7 activates the NF-kB pathway that leads to the suppression of TLR9. (A, top) C33A cells were treated with siRNA for IKKa or 
P for 16 h, and then transfected with TLR9 luciferase promoter. After 24 h, cells were exposed to the indicated treatments and harvested 24 h later to 
measure activity (left). IKKa or p levels as determined by immunoblotting (right) were analyzed. (A, bottom) C33A cells were treated with siRNA for IKKa 
or p for 24 h and TLR9 mRNA levels were measured by qPCR. Shown are data from six independent experiments performed in triplicate and error bars 
indicate SEM. (B) SiHa (HPV16+), HeLa (HPV18+), CaSki (HFV16+), and C33A control cells were treated with the IKKp inhibitor Bayl 1. At the indicated 
times, mRNA and protein expression of TLR9 as well as protein levels of kBa was determined by immunoblotting. Shown are data from seven indepen- 
dent experiments performed in triplicate and error bars indicate SEM. RE, mRNA relative expression. (C) SiHa cells were treated with Bay1 1 for the indi- 
cated time and immunofluorescence was performed to determine NF-KBp65 cellular localization (red) and TLR9 expression (green). Nuclear staining was 
controlled using DAPI. Shown are data from one out of four examined fields and one out of three independent experiments. Bars, 10 pm. (D) Human 
female primary keratinocytes (HK) were transduced with the indicated retroviruses and transfected with the NF-kB reporter gene construct. At 48 h after 
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subsequently phosphorylates and induces the ubiquitination/ 
degradation of IkBcx promoting the nuclear translocation of 
NF-kB. We noted that HPV16E7, but not E6, induced the 
phosphorylation and degradation of iKBa (unpublished data). 
To further corroborate the role of E7 on TLR9 suppression 
in model of natural HPV16 infection, SiHa cells were treated 
with a siRNA for HPV16E7 for 48 h. We observed an increase 
in TLR9 transcripts when E7 levels were suppressed. In ad- 
dition, C33A cells infected with 16QsV gained the ability to 
produce IL-8 via TLR9 stimulation only when E7 expression 
was inhibited by a siRNA (unpublished data). We then tran- 
siently expressed the NF-kB minimal promoter linked to lucif- 
erase. 24 h after transfection, cells were lysed and luminescence 
was measured. We observed that the oncoprotein E7, but not 
E6, was able to induce NF-kB activity (Fig. 2 D). Accordingly, 
HPV16 E7-transduced primary keratinocytes displayed NF-kB 
binding to a consensus cis element as shown by electromobil- 
ity shift assay (EMS A; Fig. 2 E). Antibodies against p65 and 
p50 (Fig. 2 E), but not the NF-kB family members RelB, c-rel, 
or p52 (not depicted), induced a supershift. 

We inhibited the NF-kB canonical pathway by expressing 
a nondegradable deletion mutant of IkBoi (AN-IkBo;) that lacks 
the first 36 amino acids at the N terminus containing the IKK- 
phosphorylated amino acid. AN-IkBc* expression in HPV16 E7 
HK led to cytoplasmic sequestration of NF-KBp65 and resto- 
ration ofTLR9 levels without affecting the expression of the 
viral genes (Fig. 2 F).Thus, E7 down-regulates TLR9 via the 
activation of the canonical IkBck release of NF-KBp65. 

HPV16E7 recruits an inhibitory NF-KBp50-p65 complex 
to a novel cis element on the TLR9 promoter 

The observation that NF-kB is required for TLR9 suppression 
by HPV16E7 prompted us to determine which sites of the 
TLR9 promoter are involved in this HPV16-induced event. 
Takeshita et al. (2004) previously described an NF-kB site 
that is important for the regulation of TLR9 upon its engage- 
ment by CpG 2006 (hereafter referred to as Site D). Using 
genomatix and tescan programs, we identified three additional 
NF-kB cis elements in the TLR9 promoter that we termed 
A, B, and C (Fig. 3 A). To determine which site was involved 
in TLR9 down-regulation by HPV16E7, the TLR9 pro- 
moter luciferase construct was mutated individually for the 
sites A, B, C, and D and then transiently expressed into pri- 
mary human keratinocytes transduced with HPV16E7 or the 
empty vector (pLXSN). We observed that mutation of site B, 
and not A, C, or D, restored TLR9 promoter activity in the 



presence of HPV16E7 (Fig. 3 B). Similar results were obtained 
using HK from 10 different donors transduced with HPV16 E7 
recombinant retrovirus and subsequently transfected with wild- 
type or mutant TLR9 promoter cloned in front of the lucif- 
erase reporter (unpublished data). We next determined which 
NF-kB family members bound to the site B region on the 
TLR9 promoter in the presence of HPV16E7. Chromatin im- 
munoprecipation (ChIP) for NF-KBp65, p50, PvELB, c-rel, 
and p52 was performed in HK transduced with HPV16E7 or 
stimulated with TNF or CpG 2006 (for 4 and 24 h). We ob- 
served that both NF-KBp65 and p50 forms were recruited to 
site B (Fig. 3 C) and not RelB, p52, or c-rel (not depicted). In 
addition, reChIP assays revealed that NF-kBp50 and p65 bound 
as heterodimers on site B, whereas no NF-kB complexes were 
found associated with sites A, C, or D (Fig. 3 C). In untreated 
cells an NF-KBp50 and p65 complex was constitutively re- 
cruited to site D where it most likely activates TLR9 transcrip- 
tion because high levels of TLR9 are detected in these cells 
(unpublished data). We corroborated our findings by infect- 
ing human epithelial cells with 16QsV, which also induced the 
suppressive NF-KBp50— p65 complex to site B on the TLR9 
promoter (Fig. 3 D). In addition, HPV16E7 gene silencing by 
siRNA blocked NF-KBp50-p65 binding to site B (Fig. 3 D). 
No NF-KBp50— p65 recruitment at site B was observed in cells 
treated with TNF (Fig. 3 D). 

NF-KBp50— p65 heterodimers are typically associated with 
active gene transcription (Hacker and Karin, 2006). The find- 
ing that a NF-KBp50— p65 heterodinier acts as a gene sup- 
pressor in HPV16E7-infected cells led us to examine whether 
HPV16E7 was recruited to site B along with NF-kB. To test 
this hypothesis, we transduced HK with HPV16E7HA or an 
empty vector (pBABE-HA) and performed ChIP analysis using 
anti-HA antibody. No HPV16E7-HA recruitment was ob- 
served, whereas NF-KBp50 and p65 were still found located 
at site B (Fig. 3 E). Together, these data show that HPV16E7 
induced the recruitment of an inhibitory NF-KBp50— p65 
complex at site B of TLR9 promoter without being recruited 
itself to that site. 

ERa cooperates with NF-kB complex to suppress TLR9 

As our ChIP amplification for site B is comprised of the 200 bp 
region around the cis element, we next evaluated whether the 
HPV16E7-induced inhibition of TLR9 transcription was solely 
mediated on a NF-kB cis element found on site B. We gener- 
ated an artificial promoter by cloning the site B cis element in 
front to a minimal promoter (pTAL) linked to the luciferase 



transfection, luciferase activity was determined. Data are from one assay representative of seven independent experiments; shown are mean ± SEM from 
triplicate values in three. (E) EMSA was performed on HPVl6E7-transduced HK using the NF-kB EMSA kit (Panomics) according to the manufacturer's 
instructions. For the supershift analyses, nuclear extracts were incubated with NF-KBp50 or p65 or IgG control. The arrow indicates NF-kB complexes, 
and the asterisk indicates a supershift. Shown are data from one out of four independent experiments that gave identical results. (F, left) HKs were retro- 
virally transduced with HPV16E7, empty vector (pLXSN), or HPV16E7AN-IkBo( and were harvested at the doubling population 7 and stained by immuno- 
fluorescence for TLR9 and NF-KBp65. Either Alexa Fluor 488 nm (green) or 594 nm (red) was used as secondary antibodies. Nuclear staining was 
controlled using DAPI. Shown are data from one out of six examined fields and one out of three independent experiments. Bars, 10 urn. (F, right) Western 
blot for AN-kBot (marked with an arrow), E7, and (3-actin. Shown are data from one out of five independent experiments that gave identical results. 
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Figure 3. TLR9 promoter regulation by HPV16E7 depends on the locality of the NF-kB complex at a specific cis site. (A) Predicted NF-kB 
binding sites (A, B, C, and D) in the TLR9 promoter. Sequence mutations were made at site A, B, C, and D to inactivate the binding sites. Top sequences are 
the native form; bottom sequences indicate the mutations made. (B) HKs were transduced with HPV16E7 or empty vector. Transduced cells were then 
transfected with TLR9 promoter luciferase expression vectors containing either a wt promoter sequence or sequences with the mutation in the A, B, C, 
or D site indicated in A. 40 h after transfection, luciferase activity was measured. Shown are data from six independent experiments performed in tripli- 
cate and error bars indicate SEM. (C) The binding of NF-kB complexes to the TLR9 promoter in HK was determined by ChIP and ReChIP assays. Sheared 
chromatin from HK transduced with HPV16E7 or stimulated with TNF or CpG 2006 (4 and 24 h) was immunoprecipitated with antibodies to NF-KBp50 or 
p65. Site B on the TLR9 promoter was amplified by qPCR to determine the specific binding of transcription factors bound to DNA. Immunoprecipitated 
DNA and input DNA was amplified with gene-specific and fJ-globin (Hbb) primers by qPCR, using input DNA to generate a standard curve. ChIP data are 
represented as °/o input (gene-specific)/°/o input((3-globin) = occupancy site B. Shown are data from six independent experiments performed in triplicate 
and error bars indicate SEM. (D) Re-ChIP analysis was performed for p50-p65 NF-kB complexes using C33A cells first infected with 16QsV for 36 h. Cells 
were then transfected with siRNA for HPV16E7 and harvested 24 h later. Shown are data from five independent experiments performed in triplicate and 
error bars indicate SEM. (E) E7 does not bind to the NF-kB cis element site B. ChIP was performed using human epithelial cells transfected with pbabe HA, 
pbabe HPV16E7-HA, or stimulated with CpG (4 h) or TNF (8 h). ChIP was performed using HA, NF-KBp50 or p65, or an IgG control antibody to examine 
occupancy of site B. Shown are data from six independent experiments performed in triplicate, and error bars indicate SEM. Shown is an immunoblot to 
control the ChIP for HA-E7. 
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Figure 4. The negative regulatory effect of HPV1 6 on the TLR9 promoter requires NF-kB and ERa binding to site B flanking sequences 
in human epithelial cells and HK. (A, top) Diagram illustrates the minimal B site promoter. (A, bottom) Human epithelial cervical cells were transfected 
with the NF-kB luciferase consensus site, the B minimal promoter or the control vector pTAL. 24 h after transfection 16QsV, control PV (native), CpG 
2006, or TNF were added. Cells were harvested and luciferase activity was measured after 24 h. (B, top) Diagram depicting constructs including a 200 bp 
sequence of the TLR9 promoter containing the NF-kB cis element on site B WT (B200) or mutated (Bm). (B, bottom) C33A cells were transfected with B200 or 
B200m. 24 h after transfection, 16QsV or TNF were added. Cells were harvested and luciferase activity was measured after 24 h (right). (C, left) The 200 bp 
sequence around site B on the TLR9 promoter with mapped NF-kB and ERa cis sites. (C, right) Site B was mutated at the cis ERa site (BER) and double muta- 
tion for NF-kB and ERa cis elements (BmER). (D, left) HK cells transduced with pLXSN, HPV16E7, or HPV6E7 were transfected with B200, B200m, BER, or 
BmER. 24 h after transfection luciferase activity was measured. (D, right) ChIP analyses of ERa phosphorylated and total forms and NF-KBp65 or p50 binding 
to site B on the TLR9 promoter in HK transduced ± HPV16E7. Data are representative of the mean of five or more independent experiments performed in 
triplicate, and error bars indicate SEM. **, P < 0.001 based on an unpaired Student's ttest. 



gene (termed B NF-kB; Fig. 4 A). C33A were transfected 
with the B NF-kB, control pTAL, or the NF-kB consensus 
minimal promoter linked to luciferase. Surprisingly, 16QsV 
exerted no effect on the minimal B NF-kB, although it was 
able to activate the consensus NF-kB promoter (Fig. 4 A). As 
controls, the same transfected cells were exposed to PV or TNF. 
PV did not affect the activity of any of the promoters, whereas 
TNF activated the NF-kB consensus promoter, but not site B 



NF-kB. These data suggest that additional cis elements are re- 
quired for the HPV16 inhibitory TLR9 transcriptional activity. 
To evaluate this hypothesis, we cloned a 200 bp fragment con- 
taining the site B element in the pTAL vector (B200; Fig. 4 B). 
Transfection of B200 into C33A promoted the transcription of 
the luciferase gene that was repressed in the presence of 16QsV, 
but not by TNF (Fig. 4 B). Mutations of the NF-kB site on 
B200 (Bm) restored luciferase activity (Fig. 4 B). Together, 
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Figure 5. ERa complexes with NF-KBp65 to suppress TLR9 transcription by 1 6QsV. (A, left) C33A cells were transfected with B200 and, 24 h later, 
shESRI was introduced. 12 h after shESRI introduction, cells were infected with 16QsV. Cells were harvested and luciferase activity was measured after 
24 h. (A, right) qPCR and immunoblot analysis for TLR9 and ERa in the human epithelial cervical cells treated with shESRI. (B) C33A cells were transfected 
with a shRNA scramble for ERa (left) or with an shRNAfor ERa (shESRI; right) for 24 h. Cells were then infected for 12 h with 16QsVand analyzed by 
ReChIP for NF-KBp50-ERa or NF-KBp65-ERa or NF-KBp65-p50 occupancy on site B of the TLR9 promoter. Data are representative of the mean of five or 
more independent experiments performed in triplicate; graphs show the mean ± SEM. 



these data indicate that other cis elements in the proximity 
of the predicted NF-kB cis element may be required to sup- 
press the TLR9 promoter by HPV16. We identified within 
the 200 bp region a putative estrogen response element (ERE; 
Fig. 4 C).To determine whether the ERE site was involved 
in HPV16E7-mediatedTLR9 down- regulation, we gener- 
ated mutants of B200 promoter, in which ERE cis element 
was mutagenized alone (BER) or together with the NF-kB 
cis element at site B (BmER; Fig. 4 C). Mutagenesis of ERE 
significantly alleviated the E7-induced inhibition of lucif- 
erase activity (Fig. 4 D). This phenomenon was even more 
evident in the case of BmER (Fig. 4 D, left). In agreement 
with previous data (Hasan et al., 2007a) E7 from low-risk 
HPV type 6 did not have any effect on the regulation of 
the B200, BER, or BmER promoters (Fig. 4 D, left). ChIP 
experiments in HK using antibodies against estrogen a (ERa) 
or its phosphorylated form at serine 118 confirmed that ERa 
bound to the ERE element on the TLR9 promoter in the 
presence of HPV16E7 (Fig. 4 D, right). In addition, silenc- 
ing of ERa expression by short hairpin construct against 
ERa (shESRI) restored the luciferase activity of the B200 
bp promoter in 16QsV-infected cells as well as the endog- 
enous levels of TLR 9 mRNA and protein (Fig. 5 A) . Simi- 
larly, exposure of the cells to melatonin, an inhibitor of ERa, 
blocked HPV16E7 or 16QsV ability to suppress the B200 
promoter (unpublished data). Finally, reChIP experiments 
revealed that in 16QsV-infected cells the recruitment of a 
p50— p65 complex to site B was inhibited when ERa expres- 
sion was down-regulated by shESRI (Fig. 5 B). In summary, 



our data show that ERa is recruited with the p50— p65 
NF-kB complex to site B resulting in transcriptional repres- 
sion of TLR9. 

HPV16 E6E7 induces epigenetic changes in TLR9 promoter 

Gene expression is regulated by DNA binding transcriptional 
factors and by chromatin modifications (Weng et al., 2012). 
We therefore analyzed the chromatin organization within TLR9 
promoter in mock or 16QsV-infected cells by monitoring 
Histone 4 acetylation (AceH4) and trimethylation of histone 
H3 at lysine 4 (H3K4me3), which are events associated with 
transcriptionally active chromatin (Foster et al., 2007). AceH4 
and H3K4me3 were observed at the region surrounding site 
B on TLR9 promoter in untreated C33A cells (Fig. 6 A, left). 
A similar situation was detected 6 h after infection of C33A 
cells with 16QsV, but not 8 or 36 h after infection, in which 
the AceH4 and H3K4me3 near site B were strongly decreased 
(Fig. 6 A, left). Silencing of HPV16 E7 by siRNA in 16QsV- 
C33A cells restored AceH4 and H3K4me3 association at site B 
(Fig. 6 A, right). Loss of AceH4 and H3K4me3 was not only 
limited to the TLR9 promoter region around site B but also 
occurred in the chromatin down-stream of site B until the 
transcription start site of the TLR9 promoter (Fig. 6 B). These 
histone modifications coincided with the recruitment to the 
same regions of histone deacetylases (HDACs) 1—3 (Fig. 6 C). 
Previous studies have shown that p65 can form a complex 
with HDAC3 (Xu et al, 2007). However, ChlP/reChIP ex- 
periments in cells infected with 16QsV showed that HDAC3 
recruited to the site B region was not directly associated with 
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Figure 6. 1 6QsV induces closing of the chromatin structure on the TLR9 promoter from site B until the transcription start site. (A, left) ChIP 
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the p65 (Fig. 6 D). We therefore evaluated whether ERa 
was responsible for the recruitment of HDACs to TLR9 pro- 
moter in cells infected with 16QsV. Re-ChIP experiments 



using a specific pERot (ser 118) antibody revealed that p65 
and HDAC1 were recruited in proximity to site B on the TLR9 
chromatin in 16QsV-infected C33A cells (Fig. 7 A). In addition, 
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immunoprecipation experiments with an ERa antibody 
revealed that ERa coprecipitated with NF-KBp65 and/ or 
HDAC1 in chromatin fractions extracted from 16QsV-infected 
C33A cells, whereas only a weak association of ERa in chroma- 
tin fractions from untreated C33A cells was observed (Fig. 7 B) . 
In agreement with these data, down-regulation of ERa ex- 
pression by shRNA restored AceH4 at site B on the TLR9 
promoter in 16QsV-infected cells (Fig. 7 C). Furthermore, 
inhibition of HDAC1 by Trichostatin A (TSA) restored AceH4 
and reduced HDAC1 recruitment to site B in the presence of 
16QsV (unpublished data) which coincided with an increase 
of TLR9 mRNA and protein levels. We also observed that 
16QsV infection decreases H3K4me3 at site B (Fig. 6 A). Gene 
silencing of ERa in the same cells also restored H3K4 trimethyl- 
ation near site B, indicating that ERa is also involved in 
H3K4 demethylation on the TLR9 promoter (Fig. 7 D) . We 
next wanted to determine which demethylase was involved. 
A new class of demethylase enzyme, called JARID1B, has been 
shown to be highly expressed in ERa-positive breast cancer 
cells and tumors (Dey et al., 2008; Kim et al., 2010; Catchpole 
et al, 2011; Nijwening et al, 2011). JARID1B interacts with 
ERa and catalyzes the removal of methyl groups from lysine 4 
of histone H3. We observed that JARID1B protein levels 
were elevated in HPV16E7 HK in comparison with mock- 
infected cells (Fig. 7 E, top). We hypothesized that JARID1B 
recruitment via ERa was responsible for the loss of H3K4me3 
on the TLR9 promoter at site B. Indeed, blocking ERa ex- 
pression in HPV16E7 HK decreased J ARID IB and increased 
H3K4 levels in chromatin fractions (Fig. 7 E, bottom). Re-ChIP 
experiments in 16QsV-infected C33A cells showed that 
JARID1B was recruited in association with ERa serlls at site B 
on the TLR9 promoter (Fig. 7 F). Other histone demethyl- 
ases, such as LSD1 or RBP2, were not recruited to the TLR9 
promoter (unpublished data). Abrogating ERa expression 
reduced JARID1B recruitment to site B on theTLR9 pro- 
moter (Fig. 7 G). Therefore, HPV16 induces ERa to recruit 
HDAC1 andJARIDlB histone modification enzymes as well 
as NF-KBp50— p65 to prevent TLR9 transcription. We also con- 
firmed that NF-KBp65, p50, HDAC1, and JARID1B all im- 
munoprecipitated with ERa using chromatin fractions from 
16QsV-infected C33A cells (Fig. 7 H). However, the formation 
of the different complexes were dependent on the integrity of 
the DNA because none of the subunits were immunoprecip- 
itated after DNase I treatment of the chromatin (Fig. 7 H).To 
determine whether the ERa and NF-kB complex are inde- 
pendently or dependently recruited toTLR9 promoter, we 
performed oligo pulldown experiments using biotinylated 
DNA probes which contain a region of theTLR9 promoter en- 
compassing both ERE and the NF-kB cis elements (B), intact 
ERE with a mutated NF-kB cis element (Bm), or vice versa 
(BER; Fig. 7 I). We observed that the intact site B probe se- 
quence fromTLR9 promoter (B) precipitated ERa, NF-KBp65, 
p50, HDAC1, andJARIDlB. However, mutation of NF-kB 
cis element (Bm) resulted in the loss of binding of p50 and 
p65, without significantly affecting the recruitment of ERa, 
HDAC1, andJARIDlB (Fig. 7 I). An opposite scenario was 



observed when the ERE was mutated (Fig. 7 I). Thus, the two 
repressive complexes appear to be independently recruited to 
TLR9 promoter. However, mutation of either ERE or NF-kB 
cis elements strongly affected the ERa, and NF-KBp65 inter- 
action. In summary, these data show that HPV16 promoted 
the formation of a repressive chromatin modification complex 
that negatively regulates TLR9 gene expression. 

NF-KBp65 and ERa are involved in the regulation of TLR9 
expression in HPV1 G-positive cervical cancers 

To corroborate our findings in cervical cancer samples, we ex- 
amined by irnmunohistochemical analysis the expression and 
cellular localization of NF-KBp65 and ERa in normal cervical 
tissues (n = 8) and HPV16-positive cancers (n = 8; Fig. 8). Ex- 
amination of normal cervical tissue revealed high expression of 
TLR9 in the basal (B) and suprabasal (S) layers (Fig. 8 A), which 
was lost in tumor samples (Fig. 8 B). In contrast, NF-KBp65 nu- 
clear staining was increased in tumors compared with normal 
tissue (Fig. 8, A and B). Only a subtle difference in ERa nuclear 
levels was observed in the tumor and normal tissues (Fig. 8, 
A and B). However, according to the ChIP data shown in 
Fig. 5, the NF-KBp65/ERa nuclear colocalization was signifi- 
cantly increased in tumor samples in comparison to healthy tis- 
sues (Fig. 8 C, P < 0.0001, unpaired Student's f test). To evaluate 
whether these two cellular proteins interact in cervical cancer 
cells, we performed in vivo immunoprecipitation using the 
DUOLINK technology, which determines the localization in 
tissues of two proteins in the proximity of <40 nm. With this 
technique, protein— protein interaction or proximity is revealed 
by the appearance of distinct bright dots when tissue is analyzed 
with a confocal microscope. Using specific NF-KBp65 and 
ERa antibodies, no or mainly cytoplasmic red dots (interacting 
NF-KBp65/ERa) were detected in the three normal cervical 
tissues examined (Fig. 8, D and E). In contrast, cancer specimens 
displayed high levels of NF-KBp65— ERa interactions which 
were mostly located perinuclearly or in the nucleus of the 
cells. Ortho slicer movement (total of 30 Z stack slices of 0.3 um) 
allowed us to consolidate that the red staining could be seen 
penetrating through the nucleus of the cell (Fig. 8 D, bottom). 
ChIP experiments using tissue from normal cervical tissue or 
HPV16-positive cancer tissue revealed that NF-KBp50— p65 and 
ERa were recruited to the TLR9 promoter only in cancer cells 
(Fig. 8 F). The interaction between ERa and p65 showed by 
the DUOLINK assay was also confirmed in cervical cancer cell 
lines as well as in primary keratinocytes expressing HPV16 E6 
and E7 (Fig. 8 G and not depicted). In addition, nuclear NF- 
KBp65— ERa interactions were lost in the presence of a siRNA 
for NF-KBp65 or a shRNA against ERa (Fig. 8, H and I; and 
not depicted). In summary, the analysis of human specimens 
fully confirmed the data obtained in in vitro experimental mod- 
els that showed the involvement of ERa and NF-KBp65 com- 
plex in transcriptional down-regulation of TLR9 gene. 

DISCUSSION 

The characterization of HPV mechanisms in deregulating the 
immune surveillance is extremely important to fully understand 
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Figure 7. Recruitment of epigenetic demethylating and deacetylating enzymes by ERa to site B on the TLR9 promoter in 16QsV human 
epithelial cells. (A) ReChlPfor pERa-p65 or HDAC1-3 was performed on C33A cells treated with 16QsVfor 24 h. (B, left) Immunoprecipation for ERa 
interactions with NF-KBp65 or HDAC1 was performed on chromatin fraction of C33A infected for 36 h with 16QsV. (B, right) Input controls. (C) ChIP 
using anti AceH4 histone antibodies was performed for site B on C33A cells infected with 16QsV for 24 h + shESRI. (D) ChIP using anti H3K4me3 histone 
antibodies was performed for site B on C33A cells infected with 16QsV for 24 h + shESRI. (E) Chromatin fraction Western blot analysis of JARID1 B, ERa, 
or H3K4me3 expression in pLXSN or HPV16E7 ± shESRI -transduced HK. (F) ReChlPfor pERa-p65 or pERa/HDAC1 or pERa/JARID1 B was performed on 
C33A cells treated with 16QsVfor 36 h. (G) ChIP using anti-JARID1B antibody was performed for site B on C33A cells infected with 16QsVfor 24 h ± 
shScramble control sequence (shSCR) or shESRI. (H, left) Immunoprecipation for ERa interactions with NF-KBp65, p50, JARID1 B, or HDAC1 was per- 
formed on chromatin fractions of C33A infected for 36 h with 16QsV ± DNase I treatment. (H, right) 10% of input loaded protein. (I, left) Oligo pulldown 
assay for site B, Bm, and BER using protein lysates from pLXSN or HPVl6E7-transduced cells. Bound proteins were assessed by immunoblotting for 
NF-KBp65, p50, JARID1 B, HDAC1, or ERa. (I, right) Input controls (10%). Data are representative of three or more independent experiments; graphs 
shown are the mean ± SEM from triplicate values. 
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Figure 8. HPV1 6-positive cervical cancer lesions contain NF-KBp65-ERa nuclear complexes which bind to site B on the TLR9 promoter. 

Histology and immunofluorescence (IF) of normal cervical issue (A) and HPV16+ cervical cancer biopsies (B). Nuclear (white), TLR9 (blue), ERcn (red), or 
NF-KBp65 (green). The box in the histology staining indicates which part of the slide was examined by IF. Normal (HPV— ) and a neoplastic biopsy 
(HPV16+) from one representative patient out of eight tested with similar results is shown. Bars, 10 urn. (C) Histograms representing the cellular distribu- 
tion for TLR9, ERa, and NF-KBp65 in normal (HPV—) and cancer cervical tissue (HPV16 positive by Apex screening). For each stained biopsy, six fields 
were examined and cytoplasmic or nuclear staining was counted manually and the percentage scored out of 100 cells. Data are representative of three 



JEM Vol. 210, No. 7 



1381 



JEM 



the events involved in the establishment of cervical diseases. In 
this study, we found that the oncovirus HPV16 activates a 
unique NF-KBp50— p65 and ERa inhibitory complex that sup- 
presses TLR9 transcription and function. This event resulted 
in an inhibition of IFN production, which appears to nega- 
tively influence the HPV viral life cycle. We showed that the 
oncoprotein HPV16E7 activates the NF-kB canonical pathway, 
leading to the fomiation of a suppressive NF-KBp50-p65 com- 
plex, which binds a specific NF-kB element (site B) of TLR9 
promoter. Gene silencing, chemical inhibitor, and ectopic 
mutant levels of NF-kB regulators alleviate the E7-mediated 
inhibition of TLR9 expression. Similarly, mutation of NF-kB 
site B prevented E7 to inhibit TLR9 promoter activities. Inter- 
estingly TNF, a strong NF-kB signaling activator, did not lead 
to TLR9 down-regulation. TLR9 engagement by HSV2, CpG 
oligos, or UV-inactivated 16QsV resulted in a temporary TLR9 
down-regulation, which was mediated by MyD88 and was not 
dependent on NF-KBp50— p65 binding to site B, indicating 
that one of the other identified NF-kB sites (A, C, and/ or D) 
may be required. Indeed, overexpression of MyD88DN effi- 
ciently abolished the TLR9 transcriptional repression after 
TLR9 engagement. In contrast, it did not affect the E7-TLR9 
transcriptional abrogation. These data suggest that activation of 
NF-kB signaling by different means, i.e., E7 expression, TNF 
treatment, and TLR9 engagement, led to the formation of 
distinct NF-kB complexes which may bind to exclusive sites 
in the TLR9 promoter. Furthermore, in untreated TLR9- 
expressing cells the NF-KBp50— p65 was isolated at site D 
on the TLR9 promoter, suggesting that in this context the 
NF-KBp50— p65 complex was transcriptionally active. To 
our knowledge, this is the first description of NF-KBp50— p65 
complex mediating differential regulation of a target gene 
depending on the binding site in its promoter. Generation of 
artificial minimal promoter comprising only the NF-kB cis 
element B (Fig. 4 A) revealed that additional elements were 
required to fully repress TLR9 promoter activity by E7. We 
identified, in close proximity to the NF-kB cis element at 
site B, an ERE, which we found was essential for HPV16E7 
to turn downTLR9 transcription. 



The transcription factor ERa is a member of the nuclear 
receptor family which translocates to the nucleus upon bind- 
ing to the sex hormone estradiol (Gibson and Saunders, 2012). 
Epidemiological studies showed that estrogen is a risk factor for 
both breast and HPV-mediated cervical carcinogenesis (Brake 
and Lambert, 2005; Chung et al., 2008, 2010; Shai et al., 2008; 
Chung and Lambert, 2009). In addition, experiments with 
transgenic mice expressing HPV16 E6 and E7 in the basal layers 
of the epithelia demonstrated that the estrogen cooperated with 
the viral oncoproteins in promoting cervical cancers (Chung 
et al., 2010). Our ChIP experiments confirmed that a phos- 
phorylated form of ERa (SI 1 8) was recruited to the TLR9 
promoter. However, no phosphorylation of the other amino 
acid residues known to be linked to ERa activation was ob- 
served (unpublished data). It has been previously reported that 
CDK7 is involved in phosphorylation of serll8 (Joel et al., 
1998), but it is not yet known whether this kinase is activated 
in HPV16-positive cells. In addition, blocking ERa expression 
or function in cells infected with 16QsV prevented NF-KBp50— 
p65 recruitment to site B, thus restoring TLR9 expression. Im- 
munoprecipitation of chromatin fractions from 16QsV-infected 
cells revealed the presence of ERa— p65 complexes, a complex 
which has been previously reported to act as an inhibitor for 
estrogen-regulated genes (Feldman et al., 2007). More impor- 
tantly, the ERa— p65 interaction was confirmed in HPV16- 
positive cervical cancer cell lines and tissues by DUOLINK 
assays, whereas this complex was not present in the nucleus of 
normal tissues. In addition, ChIP experiments using chroma- 
tin fractions prepared from normal cervical or cancer epithelia 
showed that ERa and p65 were recruited to TLR9 promoter 
only in cervical cancer cells. We have examined the role of ERa 
signaling on TLR9 expression in normal HK by addition of its 
ligand estradiol 17(3; indeed, we observed that TLR9 mRNA 
levels were increased (unpublished data) . These data are con- 
sistent with the nuclear expression of ERa seen in the normal 
cervical in which TLR9 expression is observed (Fig. 8). Based 
on our data, we hypothesize that ERa signaling favors cervical 
cancer development, in part by promoting an efficient and per- 
manent inhibition ofTLR9 expression only in HPV16-infected 



independent experiments performed on 16 biopsies (eight normal and eight neoplastic); graphs shown are the mean + SEM. ***, P < 0.0001. (D) Duolink 
analysis of NF-KBp65-ERa in normal versus cervical cancer tissue. DAPI marks the nucleus and the red dots represent NF-KBp65-ERa proximity ligation 
<40 nm distance. Shown are data from one out of five independent experiments that gave identical results. Represented here is one out of six fields ex- 
amined for each section. Bars, 10 um. (E) NF-KBp65 and ERa proximity interaction locality were estimated by counting the number of red dots manually 
and using the Duolink Image tool (Olink, Biosciences) present in three field sections on three normal versus three cervical cancer tissues. Data are repre- 
sentative of three independent experiments performed on six biopsies (normal and neoplastic); graphs shown are the mean ± SEM. (F) ChIP analysis was 
performed on normal cervical tissue and cervical cancer biopsies for NF-KBp65, p50, or ERa binding to site B on the TLR9 promoter (represented one experi- 
ment out of three ChIP, and qPCR for occupancy on site B was performed in duplicate three times on three normal vs. four cervical cancer tissues). Shown 
are data from six independent experiments performed in triplicate, and error bars indicate SEM. ***, P < 0.0001 based on an unpaired Student's t test. 
(G) Duolink analysis of NF-KBp65-ERa in the cervical cancer cell line CaSki. DAPI marks the nucleus and the orange dots represent NF-KBp65-ERa prox- 
imity ligation <40 nm distance. Shown are data from one out of three independent experiments that gave identical results in which one of five fields 
were examined. Bars, 20 um. (H) As in G, however cells were treated with a siRNA sequence for NF-KBp65 or small hairpin for ERa (as described in Fig. 5). 
Shown are data from one out of five independent experiments that gave identical results. Bars, 20 um. (I, left) 63x zoom x 2 and ortho Z slices (3.3 uM) show 
nuclear penetration of the NF-KBp65 and ERa proximity ligation of <40 nm in scramble short hairpin for ERa-treated cells. Shown are data from one out 
of five fields examined and one out of five independent experiments that gave identical results. Bar, 20 um. (I, right) Immunoblot demonstrating the effi- 
ciency of siRNA knockdown of NF-KBp65 and shESRL Shown are data from one out three immunoblots performed. 
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cells. Most importantly, we report that ERa was associated 
with two chromatin modification enzymes, the histone de- 
methylase JARID1B and deacetylase HDACl.This complex 
inhibited histone methylation (H3K4me3) and acetylation 
(AceH4) at site B and, consequently, downstream toward the 
transcriptional start site on the TLR9 promoter. There are an 
increasing number of reports highlighting the importance of 
deacetylases and demethylases in innate immune gene regula- 
tion. Several deacetylase enzymes, such as HDAC1, HDAC8, 
and HDAC6, influence IFN-(3 gene expression with oppos- 
ing effects (Nusinzon and Horvath, 2006). Although HDAC1 
and HDAC8 repress IFN-fi expression, HDAC6 acts as a co- 
activator essential for enhancer activity. Virus replication is en- 
hanced in HDAC6-depleted cells, demonstrating that HDAC6 is 
an essential component of innate antiviral immunity (Nusinzon 
and Horvath, 2006). We demonstrated that blocking HDAC1 
with the use of TSA restored TLR9 expression. Interestingly, 
Lin et al. (2009) treated cervical cancer cell xenografts with 
TSA and retarded tumor growth significantly. These data in- 
dicate the use of HDAC inhibitors in cervical cancer therapy 
(Takai et al., 2011). JARID1B (also known as PLU-1) has been 
shown to demethylate H3K4me3 and binds to ERa in breast 
cancer tumors (Scibetta et al., 2007; Catchpole et al., 2011). 
We show for the first time that JARID1B levels are increased 
and bind to ERa, which prevented H3K4me3 ofTLR9 in the 
chromatin fraction of HPV16-infected cells. Independently 
of ERa, JARID1B has been shown to bind to another demethyl- 
ating enzyme, LSD1, and repress the transcription of CCL14, 
an epithelial derived chemokine known to reduce the angio- 
genic and metastatic potential of breast cancer cells in vivo 
(Pedersen and Helin, 2010). These data address the role of de- 
acetylases and demethylases separately. Previous work by 
Feldman et al. (2007) showed that ERa and the p65 colocalized 
on DNA which was an essential interaction that was inhibitory 
for ER transcriptional activity. Our work further analyzed the 
ability of ERa to collectively bring deacetylase and demethyl- 
ase enzymes along with a NF-KBp50— p65 complex to silence 
the transcription of a key innate sensor. HPV16 aims to sup- 
press TLR9 as a means to avoid consequent recognition and/ 
or prevents prestimulation of TLR9 by the microflora in the 
cervix, which may be protective against HPV infection. Gillet 
et al. (2011) reported an association between alteration of the 
vaginal flora and HPV infection, suggesting that commensal 
bacterial subspecies could be protective against HPV infection 
(such as lactobacilli in the microflora) . This suppression of TLR9 
could prevent an efficient innate response against HPV and 
facilitate the establishment of a chronic infection, which is con- 
sidered a necessary condition for cervical and other virus- 
induced cancers. The importance of our findings have been 
corroborated in clinical studies showing that clearance of HPV16 
infection in women correlated to increased TLR9 expression 
in the epithelium (as well as other TLRs; Yu et al., 2010; Daud 
et al, 2011; DeCarlo et al, 2012). 

TLR9 down-regulation in HPV-induced carcinogen- 
esis is underlined by the fact that it has been shown that spe- 
cific TLR9 polymorphisms in women were associated with 



an increased risk of cervical cancer development (Roszak et al., 
2012) and that other oncoviruses, such as EBV, HCV, and 
HBV, share the property to persistendy repress TLR9 expres- 
sion although with distinct mechanism (Hasan et al., 2007a; 
Fathallah et al., 2010; Vincent et al, 2011). More specifically, 
we have recently shown that the oncoprotein LMP1 from 
EBV down-regulates the transcription of TLR9 in human 
B cells via activation of NF-kB (Fathallah et al., 2010). HBV and 
HCV particles can block the ability of pDC to produce IFN-a 
in response to TLR9, but not TLR7 agonists (Daud and Scott, 
2008; Xu et al, 2009; Hirsch et al, 2010; Daud et al., 2011; 
van Gent et al., 2011; Woltman et al., 2011). In addition, reports 
have shown that a strong dysfunction of tumor-associated 
pDCs in their capacity to produce type I IFN in response to 
CpG-A (TLR9 agonist), but not to TLR7 ligands, was observed 
in human primary breast and ovarian tumors (Hirsch et al., 
2010). These data indicate that TLR9 function is suppressed 
in viral and nonviral— associated cancers via a unique mecha- 
nism targeting TLR9 but not other TLRs which share im- 
mune signahng pathways. More work in the field of TLR9 
regulation is thus required to understand how a variety of can- 
cers affect differently the same innate immune receptor. 

In summary, our work demonstrates that the oncovirus 
HPV16 induces a transcriptional repressive complex that sup- 
presses TLR9 expression. This suggests that TLR9 may play a 
tumor-suppressive role in cervical cancers, perhaps by inducing 
type-I IFN and proinflammatory responses, which are known 
to induce cell cycle arrest, apoptosis, and death of viral infected 
cells. Thus, interfering with the regulation of TLR9 with syn- 
thetic transcriptional agonists that target ERa levels may pro- 
vide a novel therapeutic strategy for cervical cancers. 

MATERIALS AND METHODS 

Cell culture procedures. NIH3T3, Phoenix, HEK 293T, HEK 293TT 
(for virus production), and cervical cancer— derived cell lines HeLa, SiHa, 
C33A, and CaSki were maintained as previously described (Hasan et al., 
2007a). Primary human female skin keratinocytes (HK) were grown as pre- 
viously described (Hasan et al., 2007a; Mansour et al., 2007). High-titer ret- 
roviral supematants (>5 X 10 6 IU/ml) were generated as previously described 
(Hasan et al., 2007a; Mansour et al., 2007). The 16QsV and PV production, 
infection, and viral genome expression quantification of HPV16 have been 
performed as previously described (Buck et al., 2005). 

Construct information. The retroviral pBabe-puro encoding HPV16 and 6 
E6 and or E7 have been previously described (Hasan et al., 2007a). The con- 
structs pLXSN-HPV16 and HPV6 E6/E7 were a gift from D. Galloway (Fred 
Hutchinson Cancer Research Center, Seattle, WA). The plasmids used for 
HPV16 structural genes and control PV production, the target HPV16 ge- 
nome, and GFP (for PV control) were kindly donated from the laboratories of 
Martin Muller and Angel Alonso (DKFZ, Germany). The NF-kB reporter 
plasmid was obtained from BD. The TLR9 promoter luciferase construct has 
been previously described (Takeshita et al., 2004). TLR9 mutated promoters 
were generated using the Quikchange site-directed mutagenesis kit (Strata- 
gene). NF-kB minimal promoters were cloned into pTAL-LUC vector (BD). 
Minimal promoters for sites B200 and B200m were amplified from the TLR9 
promoter and mutated site B promoter, respectively, and cloned into the 
pTAL-LUC vector. The human RIG-I plasmid was donated to us from 
the laboratory of T. Taniguchi (Graduate School of Medicine and Faculty 
of Medicine, University of Tokyo, Tokyo, Japan). The AN-IkBcx, which 
lacks the sequence that codes for the first N-terminal 36 amino acids 
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(pBabe-puro-AN-lKBa), was generated by introduction of the PCR-amplified 
DNA fragment from pcDNA3-Flag-lKBa (obtained from T. Gilmore, Boston 
University, Boston, MA, and E. Kieff, Harvard Medical School, Boston, MA). 
The MyD88DN has been previously described (Hasan et al., 2004). Small hair- 
pin RNA lentiviral constructs for TLR9 (shTLR9) and control were provided 
by the Procan Axe II CLARA platform. shERa was provided by D. Picard 
(Geneva University, Geneva, Switzerland) siRNA for IKKot or IKK(3 was pur- 
chased from Ambion, and the sequence used was previously published (Accardi 
et al., 2011). siRNA for HPV16E6E7 or E7 was purchased from Dharmacon 
and previously published (Tang et al., 2006; McCloskey et al., 2010). The 
ISRE Luciferase minimal promoter was purchased from Stratagene. 

Stimuli and inhibitors. TNF (210-TA) was purchased from R&D sys- 
tems. HSV-2 was provided by the laboratory of A. Iwasaki (Yale Medical 
School, New Haven, CT). BAY 11-7082 (EMD Millipore) and melatonin 
(Sigma- Aldrich), TLR9 CpG 2006 and 2216 ODN, and negative GpC ODN 
controls, as well as RIG-I 5'ppp dsRNA were used between 3 and 10 uM 
(InvivoGen). TSA was provided by the laboratory of Z. Herceg (IARC, 
Lyon, France). The anti-IFNR and IgG control were purchased from PBL 
and used as previously described (Hasan et al., 2007b). 

Immunofluorescence. Keratinocytes were treated as previously described 
(Hasan et al., 2007a). Biopsies of normal and cancer cervical tissues were 
taken from patients, snap frozen with liquid nitrogen, and stored at — 80°C 
until required (obtained from M. Sideri, Istituto Europeo di Oncologia, Milano, 
Italia). Approval was obtained according to local ethic committees of the Istituto 



Europeo di Oncologia, Italy. Sections of 5-um thickness were cut and either 
stained for immunofluorescence using the TSA system (PerkinElmer). Cells 
or tissues were washed, the coverslips were mounted onto slides using a 1/10 
dilution of 4',6'-diamidino-2-phenylindole (nuclear stain; Invitrogen) in fluor- 
omount (Southern Biotechnology Associates), and protein expression was 
detected by direct fluorescence microscopy. Photographs were taken at 40X 
magnification. Microscopes and software used for immunofluorescence im- 
aging were performed using the NIKON Eclipse TI-NIS-Elements AR 3.10, 
Axioplan 2 epifluorescence microscope (Carl Zeiss). Data presented are rep- 
resentative images from at least four independent experiments in which >90% 
of the cells showed similar staining patterns. 

Proximity ligation duolink assay. Histological tumor sections or spotted 
cervical cancer cell lines were fixed for 30 min with cold 4% PFA. Blocking 
step was performed with Duolink Blocking solution according to manufactur- 
er's instructions. Slides were incubated overnight at 4°C with primary antibodies 
directed against NF-KBp65 (0.15 ug/ml; Cell Signaling Technology) or ERot 
(1.5 fig/ ml; Cell Signaling Technology) and then with the appropriate DNA- 
linked secondary antibodies. Duolink II Detection Orange Reagents were sub- 
sequendy used according to the manufacturer's instructions (Eurogentec). Slides 
were analyzed using an inverted confocal microscope LSM710 (Leica). 

ELISA. To measure secreted cytokines, C33A cells were seeded at 4 x 
10 4 cells/96 wells in 200 pi of growth medium. The next day the cells were 
infected with 16QsV for 24 h. Cells were washed with PBS, and CpG or 
GpC 2006 was added. 24 h later, supernatants were collected for analysis of 



Table 1. Oligo sequences used in this study 



Primers 



Forward (5'-3') 



Reverse (5'-3') 



TLR9 

P2Microglobulin 

HPV16E1 

HPV16E6 

HPV16E7 

NF-kB for ChlP 

Site A 

Site B 

Site C 

Site D 

Cloning 

Site B 

Mutations TLR9 promoter of NF-kB sites 
Site A 
Site B 
Site C 
Site D 

ERa site biotinylated 3 

NF-kB minimal promoter sites, annealing 

primers 
siRNA 

HPV16E6E7 

HPV16E7 

IKKa 

IKKfj 

NF-kB 

Scramble 



CGTCTTG AAG GCCTGGTGTTGA 
TG CTGTCTCCATGTTTG ATGTATCT 

CATAGAGATGCAGTACAGG 
CCCACAGGAGCGACCCAGAAAGTT 
As previously described (Mansour et al., 2007) 

TG G GTCTGTACCTGTGTGTG CA 
TGGATGGCCCTGTTGAGAGGG 
CTGGAGAGCACTCAGGGGAAC 
AGGCCCTGCAGAACTCTGGAG 

CCG CTAG CAG ATCTG GGGTGGGAGGTTT 

NFAF: AAGGGACTCTGGGCCCTCATCAGGCTTG 
NFBF; GAGACTTGGGGACTCGGTCAGGCAGAGGGA 
NFCF: ACA/GCG/GGT/GGA/CTT/GTC/CAT/AGG/GCC/TT 
As previously described (Fathallah et al.) 
TCAGGCAGAGGTTTCAG CACATC 
CCGCTAGCGAGTTTCTCGAGCC 



UCCAUAUGCUGUAUGUGAU 
GCACACACGUAGACAUUCG 
GCAGGCUCUUUCAGGGACA 
GGUGGAAGAGGUGGUGAGC 
As previously described (Hasan et al., 2005) 
CGAAUGUCUACGUGUGUGC 



CTGG AAG G CCTTG GTTTTAGTG A 
TCTCTG CTCCCCACCTCTAAGT 
CTCACCCCGTATAACTC 
CCCATCTCTATATACTATGCATAAATCCC 



TTCAnCCCTCCATCCACCTC 
TAGCCCCTGGGCATTCTCCTG 

GTCACACTAGGTCCCTCCTC 
TCAGGCAGAGAGCAGGGAGA 

CTCGAGCCCCTGCTTGCAGTGATCGTG 

NFAR: CAAGCCTGATGAGGGCCCAGAGTCCCTT 
NFBR: TCCCTCTGCCTGACCGAGTCCCCAAGTCTC 
NFCR: AAGGCCCTATGGACAAGTCCACCCGCTGT 

GATGTGCTGAAACCTCTGCCTGA 
G G CG ATCG CTCAAAG AGCTCGG 



a 0ther biotin-labeled oligo probes were generated using site B forward WT or mutated biotinylated primers and respective nonbiotinylated reverse primers. Probe DNA was 
amplified using the TLR9 promoter plasmid as a template. 
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IL-8, MIP3a, or IL-6 secretion using Quantikine ELISA kits (R&D Sys- 
tems) as previously described (Hasan et al., 2007a). 

Immunoblotting, immunoprecipitation, and EMSA. Biochemical 
analysis of harvested cells was performed as described previously (Hasan et al., 
2007a). To obtain cytoplasmic and nuclear extracts, cells were harvested 
and lysed as previously described (Gonda et al., 1996). Chromatin fractions 
were performed as previously described (Mendez and Stillman, 2000), omit- 
ting nuclease treatment. Where mentioned, DNase I (Fermentas) was added to 
chromatin fractions. 20 u.g of protein extracts (determined by the Bradford 
assay; Bio-Rad Laboratories) were used for immunoblotting. EMSA and 
supershift assays were performed using the NF-kB EMSA kit (Panomics). 
For each binding reaction, 5 u,g of nuclear extracts was used. Proteins or 
protein— DNA complexes were detected using ECL (GE Healthcare). Immuno- 
precipations were performed as previously described (Hasan et al., 2005). 

ChIP assay. ChIP assays were performed using the Shearing Optimization 
kit and the OneDay ChIP kit (Diagenode). For C33A cells or primary kera- 
tinocytes, cell sonication cycles last 15 s with 5 s on and 2 s off at 20% of am- 
plitude and were repeated four times. For tissue, immunoprecipitation was 
performed over night on a rotating wheel at 4°C. 2.5 ul/reaction of DNA 
solution was used for qPCR. The primers used to amplify TLR9 promoter 
regions are listed above. ReChip was performed using the diagenode proto- 
col one day ChIP kit up until step 49 and then after using the procedure 
from the Epigenome Network of Excellence website. ChIP on tissue was 
performed according to the protocol from Epigenome Network of Excel- 
lence for tissue preparation, after the Red ChIP kit from diagenode was used 
to prepare chromatin and the 1-d ChIP kit for the immunoprecipitation. 
Immunoprecipitation was performed overnight on a rotating wheel at 4°C. 
2.5 ul/reaction of DNA solution was used for qPCR. 

Chromatin fractions. Chromatin fractions were prepared as above, omit- 
ting micrococcal nuclease treatment. 

Oligo pulldown. Oligo pulldown was performed as previously described 
(Lopez-Rovira et al., 2002) with nuclear extracts as stated in the figure leg- 
end and oligo probes as listed in Table 1. 

Transfections and luciferase assay. Cells were transiently transfected 
with the luciferase constructs or sh vectors using FuGene (Roche) as de- 
scribed previously (Hasan et al., 2007a). Each experiment was repeated three 
times in triplicate; results generally deviated by <10% of the mean value. 
SiRNA were transfected as previously described (Hasan et al., 2005). 

Type I IFN bioassay. Supernatants were harvested, UV inactivated, and 
placed onto transfected HEK293T cells that express the IFN-fS— inducible cis 
element ISRE-linked to the luciferase gene. 24 h after stimulation with su- 
pernatants, cells were harvested and luciferase activity was measured as previ- 
ously described (Hasan et al., 2005). 

Genotyping. Tumor samples were genotyped using multiplex PCR with 
HPV type-specific primers for amplification of viral DNA and array primer 
extension for typing (Hasan et al., 2007a). 

RT-qPCR. Total RNA was extracted from cells using the RNeasy Mini kit 
(QIAGEN and Machery Nagel). cDNA was synthesized with the First strand 
cDNA synthesis kit (MBI, Fermentas). The Mx3000P real-time PCR sys- 
tem (Stratagene) was used to perform qPCR with Mesa green qPCR Master- 
Mix plus (Eurogentec) . Primer sequences are enclosed in Table 1. 

Statistical analysis. GraphPad (version 5) was used to calculate unpaired 
and paired p-values. 
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